Spin relaxation in a Rashba semiconductor in an electric field 
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The impact of an external electric field on the spin relaxation in a disordered, two-dimensional 
electron system is studied within the framework of a field-theoretical formulation. Generalized 
Bloch-equations for the diffusion and the decay of an initial magnetization are obtained. The 
equations are applied to the investigation of spin relaxation processes in an electric field. 
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I. INTRODUCTION 

In recent years much attention has been paid to the 
emerging field of spintronics in semiconductor physics. 
Much of the interest in this field is stimulated by the 
spin-ficld-cffcct transistor proposed by Datta and Dasi. 
In this proposal the Rashba spin-orbit interaction is used 
to manipulate the charge transfer through a transistor 
like device by controlling the electron spin. The Rashba 
interaction, which is the conventional spin-orbit interac- 
tion with a constant electric field, affects the electron spin 
like a momentum dependent magnetic field. 

In practice the realization of the proposal by Datta and 
Das has proven to be a challenging task. In particular 
the injection of spins into non-magnetic semiconductors 
has turned out to be demanding, although successful in- 
jection has been reported by several authors^ 3 * 4 *^. De- 
spite this fact the importance of the Rashba interaction 
for the concept of the spin-field effect transistor has mo- 
tivated further studies on its impact on charge and spin 
transport properties. Particular attention has been paid, 
e.g., to investigations of spin-accumulations 7 - and to the 
spin-galvanic effect, to an intrinsic spin-Hall effect, and 
to the weak localization-antilocalization transition in the 
conductivity 10 . 11 ! 12 . 13 ! 14 . 

Since in the Rashba model the electron wave-vector 
is coupled to the spin the spin is quickly randomized 
by scattering events. Therefore, the ballistic transport 
regime is the most important regime for the spin-field 
effect transistor proposed by Datta and Das. Accord- 
ingly, most of the investigations focus on the situation 
in which the disorder energy is small compared to the 
Rashba level splitting. Recently, however, a generalized 
spin-field effect transistor has been proposed^ which is 
also expected to work in the diffusive regime, even if the 
Rashba level splitting is small compared to the charac- 
teristic disorder energy. This raises the question which 
impact of the Rashba interaction on the charge and spin 
transport properties can be expected in this limit. 

The aim of the paper is to provide a theory for the spin- 
relaxation in a a two-dimensional electron gas in the dif- 
fusive regime, in the limit that the characteristic disorder 



energy is large compared to the Rashba level splitting on 
a time scale, in which inelastic processes can be ignored. 
To this end we generalize the non-nonlinear cr-model to 
include also the impact of the Rashba interaction and 
an in-plane electric field. Using the generalized cr-model 
we derive generalized Bloch-equation and use these equa- 
tions for the investigation of an initial magnetization. 



II. BASIC EQUATIONS 

We consider charge carriers in a two-dimensional plane. 
The motion of the charge carriers is described by the 
Hamilton operator 

P 2 

H = — + F- x + a ■ {N x p) +V(x). (2.1) 
2m 

Here p is the momentum operator, m is the effective 
mass, F is a constant field in the two-dimensional plane, 
AT is a constant vector perpendicular to the plane, and 
V(x) is a random potential. The components of the vec- 
tor er = u x e x + cr y e v + a z e z are the Pauli matrices. e x 
and e y are the unit vectors in the two-dimensional plane 
and e z is a unit vector transverse to the plane. Accord- 
ingly, N = N z e z and F = Fe x . The random potential 
is characterized by a Gaussian distribution function with 
zero mean and second moment 



(V{x)V(x')) di 



1 



■5{x-x'). 



(2.2) 
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Here (. . . )dis denotes the configuration average, v is the 
density of states per spin at the Fermi- level for N z — 0, 
and r is the single-particle scattering time. 

In our investigation we focus on configuration averages 
of retarded and advanced Green functions. These func- 
tions are solutions to the differential equation 

((±ioj + E)l - H)G r/a (x,x'\E,lu) = 16(x-x'). (2.3) 

In this equation G R I A is a 2 x 2 matrix in spin-space , 1 is 
the 2x2 unit matrix and E and uo are fixed energies. As a 
boundary condition we require that the Green functions 
vanish at infinity. In this case the configuration averaged 
Green functions are invariant against translations by a 
fixed vector a if the energy E is also shifted by F ■ a. 
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Accordingly, they satisfy the relationship 

( x + a, x' + a\E, u>) 
= (G R / A ) dis (x,x'\E-F-a,uj). (2.4) 

To calculate the configuration averages of the Green func- 
tionswe use the field theoretical formulation of the Refs. 
|!6ll7ll8| . In this formulation the replica trick is used 
to derive an effective action for the calculation of the 
configuration averaged Green functions. Following the 
standard derivation we find that for the problem at hand 
the partition function takes the form 

Z = J DQex.p(A Q ), (2.5) 



where 



Ac 



-^trQ^trlnG^ (2.6) 



Here Q is a hermitian (2n r + 4) x (2n r + 4) matrix field 
(n r -numbcr of replica copies+ 2 components particle- 
hole space-t- 2 components spin indices) with compo- 
nents Q^ m a i n { x ) (a, a'- replica indices, i, j- particle-hole 
indices, m,n-spin indices). It satisfies the relationship 
QC = CQ T , where Q T is the transposed matrix and 



C = 



l<7y 

ia„ 



(2.7) 



The tr in Eq. H2.6[) contains both the summation over the 
discrete indices and the integration over space. Further- 
more, 



G 



Q 



(x, x') = [ {(iiv a + E)S mn ~ H mn }S aa 'Si 3 



III. THE SADDLE-POINT GREEN FUNCTIONS 



To calculate the configuration-averaged Green function 
in saddle-point approximation we now look on the saddle- 
point solution of the effective action Aq , defined by 



SAq 
6Q 



= 0. 



(3.1) 



The saddle-point values of Q and Gq are diagonal with 
respect to the replica indices, as are the exact expec- 
tion values. Accordingly QsPaLa'n = Qamn^aa' and 
G Q P Z a >n = G «iJ aa -- From the Eqs. (H and O 
we find 



Q Qr l(x) = — G a % n (0,0\E-F-x,Lj) 



Here fi x = E — F ■ x. 



(3.2) 



To find an equation for the Green function we use 
the gradient expansion 1 -. To this end we first introduce 
Wigner coordinates, according to R = (x + x')/2 and 
r = x — x' , and perform a Fourier transformation with 
respect to the relative coordinates. The Green function 
in the new coordinates G is related to the old Green func- 
tion by the relationship 



G a Z n (R + r/2,R-r/2\E,u) 



dk 



e* k - r G a % n (R,k\E,iu). (3.3) 



where co a — u;A a , A a = 1 for a > and A Q = — 1 for 
a < 0. The index a can take on values between — (n r — 1) 
and n r . It is understood that the limit n r — > is taken 
at the end of the calculation. 



In calculating G we restrict the consideration to the low- 
est order in the gradient expansion. In this approxima- 
tion the equation for the Green function takes the simple 
form 



J 



j 



{[(iu a + hr - e fc )l — ha ■ (TV x fc)]<% + i—Q 1 £^r))G£(R, k\E,w) = IS, 



2t 



(3.4) 



r 



Here all objects are considered as 2 x 2 matrices in spin 
space, e fc = h 2 k 2 /2m and fi R = (i x \ x =r- 

In order to solve the Ecis. H3.2| l- H3.4f> self-consistently we 
restrict the consideration to the classical accessible region 
fj,R > 0. Furthermore, we consider only the dirty limit, 
in which both the disorder energy h/2r and /ir are large 
compared to the Rashba level splitting \N z \x/2miJ J R and 



assume that \xr is large compared to h/2r. In this case 
we find that the self-consistent solution to the equations 
(13.21) and (|3.4|) takes the form 



6 mn SijSgiiu! a , 



(3.5) 
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G a % n {R,k\E,u>) 



(k\/j, R ,oj)Sij, 



(3.6) 



g a (k\fiR,u) = 



M a (k\^ R ,w) 



J a + f-lR - £k + l-gpSgnUJ a 



iN,Hk 



+ 



l0J o 



-iN z hk- 
- f-k + ' 



£sgnw 



(3.7) 



M a (k\n R , lu) 



[i(oj a + — sgnw Q 

ZT 



Hr 



E + (k))(i(uj a + — sgnuj a ) + fi R 
It 



E-(k)), 



(3.8) 



r 



where k± = k x ± ik y and E± — ± h\N z \k. In writ- 
ing down the saddle-point field 13.5fl we have ignored 
its imaginary part which only leads to an a-independent 
shift of hr. We would like to note that the simple struc- 
ture of the saddle-point field is a consequence of the fact 
that the density of states for N z = is constant. In an en- 
ergy dependent density of states the saddle-point would 
pick up an energy and field dependent contribution from 
the density of states, which would modify Ea. (|3.5|) . 



Such matrices are generated by hcrmitian generators W 
(B = exp(iW)), which have the structure 



A=0 i=l 

+S aal lcT3®Ta + iA aa , 3 l a 3 ®(7 i (4.2) 



IV. THE NON-LINEAR cr-MODEL 

The fluctuations around the saddle-point solution are 
characterized by massive and soft modes. The latter de- 
scribe the diffusion and spin relaxation properties we are 
interested in. In the field-theoretical formulation the ex- 
istence of these modes is a consequence of a continuous 
symmetry of the effective action (|2.6fl which exists in the 
limit uj, N z — * 0. In this limit the effective action is in- 
variant against similarity transformations with hermitian 
matrices B, which satisfy the relationship 



B 1 CB = C. 



(4.1) 



Here S aa '^ (A aa >^) are real matrices symmetric (anti- 
symmetric) in a and a', oo = 1, and i = 1,2,3 corre- 
sponds to x, y and z, respectively. The components of 
the matrices W aa ' , which are associated with the upper 
index of the matrices S aa >^ and A aa /^ (that is with A), 
correspond to the particle-hole indices. The lower indices 
characterize the spin. 

In order to derive an effective action for the soft 
modes we follow the standard derivation for the non- 
linear sigma-modeli£iiL 1 £. The complications originat- 
ing from the presence of the electric field we treat as in 
Ref . [20| . Doing so, we find that the effective action for 
the diffusion and spin relaxation modes takes the form 



dx9(fi x )tr(AQ{x)) 



dx6([i x )D(ii x )tr(DQ(x) ■ DQ(x)). 



(4.3) 



r 



Here Q(x) is the soft part of the matrix Q(x) generated 
from the saddle-point field Qsp(x) by the local similarity 
transformation Q(x) = B(x)Qsp(x)B~ 1 (x), 



(4.4) 



The bracket [...,...] symbolizes the commutator and 
the step function 0(fi x ) restricts the range of integration 
to the classical accessible region. The matrix er X N in 
Eq. (|4.5() is a matrix which only acts on the spin- indices 
but not on the replica and the particle- hole indices. 



and 



DQ(x) = VQ(x)-i—[ax N,Q(x)}. (4.5) 



In order to investigate the non-linear cr-model further 
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we write the matrix Q in the form 

Q=( ^ (4.6) 

V 921 ~~ V 1 - 521512 / 

where 

<7i2 = q aa >\V K ® cta, (4.7) 

a > 0, a' < and 921 = 9i2' wnere + symbolizes her- 
mitian conjugation, and expand the action i|4.3[) with re- 
spect to powers of q up to fourth order in q. Doing so, 
we use the fact that the matrices q a a'ii are antisymmet- 
ric with respect to a and a' for A = 0, 1, 2 and n = and 
for A = 3 and k — 1,2,3 and symmetric with respect to 
a and a' for A = 3 and k = and for A = 0, 1,2 and 
k = 1,2,3. The quadratic terms in the expansion yield 
the Gaussian propagator, the physics of which agrees 
with the ladder-approximation. The quartic terms in the 
expansion yield the weak-localization corrections to the 
transport coefficients. 



^s{fJ-x) = —D(fj, x ), (4.13) 

and Di(/j, x ) = D(p x ). 

According to the Eas. (|4.9|) . I4.1L)[I and 14.1111 the equa- 
tion for the particle transport decouples from the equa- 
tions for the spin transport. The particle transport is 
described by the propagator Poo- Its equation of motion 
contains besides ordinary particle diffusion also the par- 
ticle packet drift and the heating of the charge carriers. 
Its physics is detailed in Ref . |20j| . The spatial compo- 
nents of the diffusion propagator, for i,k — 1,2,3, 
describe spin diffusion and relaxation. In addition to 
the diffusion, drift and heating processes present in the 
particle number diffusion equation, these equations also 
contain the spin relaxation and an electric-field depen- 
dent rotation of the spin-direction, which is induced by 
the Rashba-interaction. We would like to note that the 
same equations can also be obtained from the ladder ap- 
proximation. 



A. Diffusion and relaxation in the Gaussian 
approximation 

In the Gaussian approximation the effective action 
takes the form 



»>° A.k,<5" 

'<o ' 



x q aa '1(x)T XS {x\E,s)ri^q aa ^{x),(A.8) 

where 77^ = 1 for A = 0, 1, 2 and k = 1,2,3 and for A = 3 
and k = and = —1 else. The operator F\g(x\E, s), 
the non-vanishing components of which are 



T 00 (x\E, S ) = s-(V,D(^)V) (4.9) 



and 



r ik (x\E,s) = {s + n i (n x )(l + 5 i3 )-{V > D i (n x )V))5 ik 
+ -^-{Sk3{^»{fix), Vi} - {uj s (n x ),Vi}8 i3 ), 

(4.10) 

is related to the generalized diffusion propagator 
P\ K (x,x'\E,s) by the relationship 

r KX {x\E, s)P X5 {x, x'\E, s) = S kS S(x - x 1 ). (4.11) 



Here the bracket {.., ..} symbolizes the anti-commutator, 
hs = lo/2 and the indices i and k can take on the values 
1, 2, 3, which correspond to x, y, and z, respectively. The 
transport coefficients are given by 



4N?m 2 



r, 2 



--D(Mx), 



(4.12) 



B. Quantum corrections 

In order to calculate the quantum corrections to the 
transport coefficients we also take into account the quar- 
tic terms in the expansion of A a with respect to q and 
investigate the corresponding functional integral in one- 
loop approximation. Doing so, we restrict the consid- 
eration to the linear response regime. In this case we 
can ignore the impact of the electric field on the trans- 
port coefficients and thus set the external electric field F 
equal to zero in all loops. For charge transport the range 
of validity of this approximation has been discussed in 
Ref. [2(J • There it has been shown that the impact of the 
field is negligible, as long as the spread of the particle 
packet is small compared to the distance of its center 
from the turning point. The time needed by a particle 
packet, which was a delta-pulse centered at x = at time 
t = and in which all particles have the same energy E, 
to reach a width of the order of E/F is for a typical ex- 
periment larger than 10~ 4 s, as shown in Ref. |20l| . Given 
that the typical lifetime of an initial magnetization is in 
the nanosecond range it seems to be reasonable to restrict 
the consideration to this limit. 

To one-loop order the particle number diffusion coeffi- 
cient D(E) in Ea. H4.9(l is replaced by 

1 



dW(E,s) = D{E)(1 



2-Kvh 



(P 33 (E,s)+2P 11 {E,s) 
-P 00 (E,s))), (4.14) 



where 



Poo(E,s) 



1 

4ttD(E) 



In 



D(E)A 2 



P n (E,s) 



4irD(E) 



In 



D{E)A 2 
s + Sl(EY 



(4.15) 



(4.16) 
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and 



^33 (E,S) 



1 



In- 



D(E)A 2 



A-kD{E) ^ s 2 + 29, 2 {E) +3sn(E)' 

(4.17) 

Here A is the ultraviolet cutoff. These corrections have 
the same structure as the conventional anti-weak local- 
ization corrections to the conductivity^. In a current 
relaxation experiment the frequency s is effectively re- 
placed by l/r^(T), where r^(T) is the inelastic phase 
breaking time which increases with decreasing tempera- 
ture. Accordingly, the particle number diffusion coeffi- 
cient has a minimum as a function of the temperature at 
T = T , where fl(E) m I/t^Tq). Therefore, the diffu- 
sion coefficient increases with decreasing temperature if 
T < Tq. 

The opposite turns out to be true for the spin-transport 
coefficients. To one-loop order the spin diffusion coeffi- 
cients D X (E) and D y (E) are replaced by 



D^ v (E,s)=D(E)(l 
D Z {E) by 



1 



2-Kvh 



(P 00 (E,s)+P 33 (E,s))), 
(4.18) 



dP(E,s)=D(E)(1--^—(P 00 (E,s) + 2P 11 (E, S ) 
Zirvn 

-P 33 (E,s))), (4.19) 



the spin relaxation frequencies Qi(E) by 



n x/y {E,s) = tl(E) 



D { P{E,s) 
D(E) 



and 



n^(E,s) = n(E) 



D X 1 \E, s) 
D(E) 



(4.20) 



(4.21) 



and the rotation frequency uj s (E) by 



wW(£,«)=u.(£)(l 



1 



2nvh 



(P 00 (E,s)+P n (E,s))). 

(4.22) 

In these equations s has also to be replaced by 1/t^(T) 
if s is smaller than 1/t^T). In contrast to the particle- 
number diffusion coefficient, however, the spin-transport 
coefficients have no minimum as a function of the tem- 
perature. In the range of applicability of the Eas. j4.18jl - 
l|4~2*2*|l the magnitude of the spin transport coefficients 
keeps decreasing with decreasing temperature even if 
T < Tq. In contrast to the corrections to the particle- 
number diffusion coefficient, which are anti- localizing at 
low temperatures, the corrections to the spin transport 
coefficients always show a tendency to localization. 

At this point we would also like to mention that 
the quantum corrections render the spin-transport co- 
efficients weakly anisotropic. This anisotropy results 



from the fact that Pn(E,s) ^ P 33 (E,s). However, 
since the difference between these quantities is of the or- 
der of l/A-KvhD(E) this difference is small compared to 
the magnitude of the quantum corrections themselves, 
which are larger than the former by a factor of the or- 
der of \n{D(E)K 2 / (s + Q(E))). Therefore, we ignore 
the anisotropy in the following and replace P33 (E, s) 

by Pn(E,s) everywhere. In this case Cl^iE, s) — 
S$\E,s) = flW (£?,«) and D x 1] (E,s) = D^(E,s) = 
D^s\e, s) 



V. APPLICATION TO RELAXATION 
PHENOMENA 

The evolution of an initial spin packet is described by 
the Bloch equations. Phenomenology tells us that the 
Bloch equations have the structure 

sS(s) + n(s)- S(s) + R(s)x S(s) + divJ(s) = S . (5.1) 

Here S(s) is the Laplace transform of the spin density 
with respect to time, s is the corresponding Laplace 
frequency, Sq is the initial condition, fi(s) is a second 
rank tensor, R(s) is a vector, and J(s) is the spin- 
current tensor. For the components we use the nota- 
tion (f2 • S)i = J2k ^ikSk and (divJ)i = J2k ^kJki- The 
quantities S(s), So and J(s) depend also on x and on the 
total energy E, but this dependence has been suppressed 
in En. lj5.ip to simplify the notation. Due to the general- 
ized translation invariance (j2.4(l the transport coefficients 
i7 and R do not depend on the quantities x and E sep- 
arately, but only on the kinetic energy [i x . In Ea. (|5.1() 
this dependence has also been suppressed to simplify the 
notation. 

From the physical point of view the tensor Sl(s) de- 
scribes the decay of the initial magnetization and R(s) 
describes the rotation of the initial magnetization. The 
disorder renders both quantities dispersive. Therefore, 
both n(s) and R(s) depend on s. This dependence 
should be observable in optical experiments where it 
should reflect itself in deviations from the purely expo- 
nential decay. 

The equations (j4.10|l and Ij4.11jl can be cast into the 
form Ij5.1jl . Doing so, it has to be taken into account 
that the bare transport coefficients are renormalized by 
the one-loop corrections of the preceding section. From 
Eci. (j4.10JI we find 



and 



n«(Mx,«) = n (1) (p x ,s)S tk (l + 5 a ), (5.2) 



(5-3) 



2 dfi x 

The spin current tensor takes the form 
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J ki {x\E, s) = -D^{fi x , s)V k Si{x\E, s) + uj^\ f i x ,s){N l S k (x\E, s) - 5 lk (N ■ S(x\E, a))). 

I 



(5.4) 



This tensor has been obtained by collecting all those 
terms in Ea. H4.10fl . which can be cast into a divergence. 
Accordingly, the tensor l|5.4|l might deviate from the true 
spin-current tensor by terms which are annihilated by the 
divergence. However, since every observable quantity can 
be calculated directly from Ea. (|5.1|l such terms are not 
important for our further considerations. 

The Eas- H4.14fl - H4.22fl ignore the impact of inelastic 
scattering events on the transport coefficients. To take 
them into account we introduce the phenomenological 
phase relaxation time (T) and replace s — > s + 1 /r^ (T) . 
In this case the transport coefficients are given by the 
Eas. lfS~2|) and (|5"3|) with s replaced by 1/t^(T) for s < 
1/t0(T). Accordingly, the relaxation becomes Markovian 
for t > T0 (T) . Below we restrict the consideration to this 
situation. 



A. Decay of a homogeneous initial magnetization 



0.4 
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In order to investigate some of the consequences of 
Ea. (|5.1|) we first focus on the decay of a spatially homo- 
geneous magnetization. To this end we assume that an 
initial magnetization So has been created at time t = 
and that only charge carriers in the vicinity of the Fermi- 
energy fi contribute to the magnetization. In this case the 
initial distribution function is given by 



S°\ F =o(x,E) = S 6(fi-E), 



(5.5) 



in the absence of the field. In the presence of the field 
the energy of a charge carrier placed at x is changed by 
F ■ x. Therefore, the initial condition takes the form 



S°{x,E) = S S(ji + F ■ x - E) 
= S S(fi- fi x ) 



(5.6) 



in the presence of the field. 

During the evolution the spin distribution remains spa- 
tially homogeneous. Accordingly, S(x\E,s) = S(fi x ,s). 
The direction and the magnitude of the magnetization 
and the energy of the particles are changed. Since, how- 
ever, the magnetization decays on a time scale small com- 
pared to the time needed by a particle to change its en- 
ergy appreciably the energy distribution function remains 
close to the initial delta-function during the time in which 
the magnetization is measurable. To take advantage of 
this fact we make the ansatz 



S(fi x , s) = S(s)8{fi - fi x ) 



and derive a closed equation for S(s) by integrating Eq. 
(|5.1|) with respect to fi x . Doing so, we obtain the simple 



FIG. 1: Decay of S z for 2fi F /^ (1) = 0.3 (solid line), 2 (dashed 
line) and 10 (dashed dotted line). The initial conditions are 
S z (0) = 1 and S x (0) = 0. 



algebraic system of equations 

(,s+^ 1 )( M ))^( S )-i^0^[(iVxF)x5( S )] ;E/ , = S° x/y , 

(5.8) 



(s + 2rt 1 H f i))S z ( S ) l^^l( N >< F ) >< = S° z . 

(5.9) 

This system has the same structure as that for the decay 
of an initial magnetization in an ordered system (hop- 
ping transport of small polarons on an ordered lattice)^. 
Only the structure of the transport coefficients is differ- 
ent. 

To calculate the magnetization we solve this system 
and perform an inverse Laplace transformation. Doing 
so, we obtain 



S y (t) = S°e- fi(1 \ 



(5.10) 



(5.7) ^(t)= e - 3 ° (1)t / 2 [ S° x [co S h(cu F t) + 



fiW sinh(wFf) 



o sinh^F*)-, 



(5.11) 
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FIG. 2: S x (t) for S" 2 (0) = 1 and 5,(0) = 0. The parameters 
are the same as in Fig.l. 



and 

a -3£2 (1) t/2 r oOr w ^ O^ 1 ' sinh(wp£) , 

S z (t) = e iU t/z [ S z [cosh.(ui F t) -] 



UJf 



Here O^ = O^^), 



UJF 



4J2I./2, 



and 



2 



(5.12) 



(5.13) 



(5.14) 



The Eas. (|5~TUjl - l(5~T2ll apply for t > r^(T). They do not 
contain the magneto-electric effect of Ref. ^23] , since this 
effect is proportional to the ratio between the Rashba 
level splitting and the Fermi energy, which is considered 
as small in our treatment. 

According to the Eas. (|5.1U|l - H5.14|l the character of the 
decay of the initial magnetization depends on the mag- 
nitude of the electric field. If F < F c , where F c is deter- 
mined from the requirement O^ = 2|fiF c |, the frequency 
ujf is r eal and reflects itself in an exponential decay of 
an initial magnetization (see the Figs.l and 2). Further- 
more, the i-component of the magnetization is coupled 
with the z-component of the magnetization. If S z ^ 
a non vanishing magnetization in x-dircction is created 
even if 5° = 0. If F > F c the frequency u>f becomes 
purely imaginary. In this case the hyperbolic functions 
turn to trigonometric functions, which lead to a rotation 
of the decaying initial magnetization (see Fig.l and 2). 



Therefore, the vector (S x (t), S z (t)), with components 

fiW sm(\uj F \t), 



S x (t) 



_ g -3n (1) t/2 



S2[C0B(|W*-It) 



\UJ F \ 



sin(|u/i? |t). 



(5.15) 



and 



S z (t)=e- 3n(lH / 2 [ S°[cos(K^) 



ftW sm(\u) F \t), 



UJp 



sin(|o; F |t) 
-b x il F — r— j — j 
\uj f \ 



(5.16) 



moves on an ellipse with decaying diameter. The ellipse 
satisfies the equation 

n F S x 2 (t) + nS x (t)S z (t) + Q F S 2 z (t) = const-exp(-30 (1) f) 

(5-17) 

The quadratic form can be diagonalized. Doing so, it can 
be shown that the elliptic axis are tilted by 7r/4. 



B. Alternating Hall-current 

The rotation of the magnetization discussed in the pre- 
ceding subsection reflects also in the current. In the pres- 
ence of the Rashba interaction the current operator has 
the form 2 ^2£ 



p{t) 2e 



S(t) x TV. 



(5.18) 



Here e is the charge and S is the spin operator. The 
Heisenberg equation for p(t) is independent of S(t). 
Accordingly, the configuration average of the expection 
value of the first term on the right hand side of Eq. I|5.18|) 
yields the conventional Drude current and its weak lo- 
calization corrections. The configuration average of the 
second term yields two additional contributions, a cor- 
rection term to the longitudinal current 



Sjx(t) = ~N z Sy{t), 



and the Hall-like term 
hit) 



2c 



N z S z (t) 



(5.19) 



(5.20) 



perpendicular to both F and N. If we chose the same 
initial condition for the magnetization as in the preceding 
section S z (t) is oscillating, if 2\Op\ > OS 1 '. According to 
Ea. (|5.2()fl these oscillations manifest themselves in an al- 
ternating Hall-like current. The alternating Hall current 
is produced by a constant electric field. However, since 
the magnetization decays when time elapses the current 
oscillations are not stable but decay too. This raises the 
interesting question whether these oscillations can be sta- 
bilized by some means. 
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VI. RESULTS 

In the paper we have used field theoretical methods 
to derive the generalized Bloch equations for the mag- 
netization in a disordered two-dimensional electron gas. 
Our equations apply for time scales which are large com- 
pared to the elastic scattering time, large compared to 
the phase-breaking time but small compared to energy 
relaxation time. Accordingly, they take into account the 
particle- and spin diffusion, the drift induced by an ex- 
ternal electric field in the plane, the heating of the charge 
carriers, the decay of the magnetization due to scatter- 
ing events, and an electric field dependent rotation of the 
magnetization direction. 

The calculation of the transport coefficients shows that 
the impact of the quantum-corrections on the spin trans- 
port coefficients is different from that on the particle 
transport coefficients. The corrections to the latter are 
antilocalizing. This fact manifests itself in a minimum of 
the diffusion coefficient as a function of the temperature. 
The structure of the corrections to the particle trans- 
port coefficients agrees with the conventional antilocal- 
ization corrections to the conductivity due to spin-orbit 
coupling 21 . In contrast to the particle transport coeffi- 
cients the spin transport coefficients show a tendency to 
localization. Both the spin diffusion coefficient and the 
spin decay rate decrease with decreasing temperature, 
even at very low temperatures. This sets the situation 
for the spin-transport coefficients apart from that for the 
particle transport coefficients, which increase with de- 
creasing temperature at sufficiently low temperature due 
to antilocalization. 

We have applied these equations to an investigation of 
the decay of a spatially homogeneous initial magnetiza- 
tion. In the absence of the in-plane electric field such a 
magnetization decays exponentially at large tiniest. The 



decay is anisotropic. Its transverse component decays 
twice as fast as its longitudinal components. The temper- 
ature dependence of the decay is given by the tempera- 
ture dependence of the weak localization corrections to 
the transport coefficients, which lead to a logarithmic re- 
duction of the decay rate with decreasing temperature. 
Accordingly, the decay becomes slower with decreasing 
temperature. 

Already a weak electric field renders the decay qualita- 
tively different. If the electric field is smaller than a crit- 
ical field it tips the magnetization only slightly over. If, 
however, the field exceeds the critical field F c the magne- 
tization starts rotating in the plane transverse to TV x F . 
The rotation frequency depends on the degree of disor- 
der, on the magnitude of the Rashba-parameter, on the 
magnitude of the electric field and on temperature. An 
increase of the temperature results also in an increase of 
the rotation frequency lof and in an increase of the criti- 
cal field F c . We would like to note that a similar rotation 
has also been obtained in Ref.[2^] for small polarons on 
an ordered lattice. Our calculation shows that this ro- 
tation can also be observed in a disordered electron gas, 
in which the transport mechanism is very different from 
that for a small polaron on an ordered lattice. 
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